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ABSTRACT 



The phase response of a network is measured at uniform 
frequency intervals. A linear regression analysis is per- 
formed on samples of the phase response measured at 
frequencies within an aperture centered on a group delay 
frequency to obtain an estimate of the group delay of the 
network at that frequency. The process is repeated for a 
sequence of group delay frequencies to determine a trace of 
the group delay of the network across a range of frequencies: 

20 Claims, 4 Drawing Sheets 
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GROUP DELAY ESTIMATE SYSTEM USING 
LEAST SQUARE FIT TO PHASE RESPONSE 
RAMP 

this is a divisional of application Ser. No. 08/013,882 
filed on Feb. 5, 1993, now U.S. Pat. No. 5,399,976. 

FIELD OF THE INVENTION 

The present invention relates to electronic measurement 
instruments, and more particularly to a method and appara- 
tus for estimating the group delay of a network with such an 
instrument. 

BACKGROUND AND SUMMARY OF THE 
INVENTION 

For expository convenience, the present invention will be 
illustrated with reference to one particular application 
thereof, namely as a group delay estimator in a network 
analyzer. However, it should be recognized that the inven- 
tion is not so limited. 

In the field of the electrical arts, it is often necessary to 
analyze the electrical characteristics of a configuration of 
interconnected electrical components or "network." One 
characteristic of networks that is often of interest is the 
length of time required for a signal to travel from an input 
to an output. A measurement made on networks and related 
to this characteristic is known as the group delay. The group 
delay of a network is defined as the negative of the deriva- 
tive of the phase response of the network with respect to 
frequency or, in other words, the rate at which the phase 
response changes with frequency. 

Traditionally, the group delay of a network has been 
determined by measuring the phase response of the network 
at two different frequencies using a stimulus signal applied 
to an input of the network. The difference between the two 
frequencies is termed the "aperture." The group delay of the 
network at a frequency in the center of the aperture is 
calculated as the difference in unwrapped or cumulative 
phase between the two phase response measurements 
divided by the aperture width. Usually, the phase difference 
is measured in radians and the aperture width is measured in 
radians per second to yield a group delay measurement in 
seconds. However, other units can be used. 

Previously, group delay measurements required special 
purpose, dedicated hardware. More recently however, elec- 
tronic measurement instruments known as network analyz- 
ers have been used. Network analyzers can be used to 
provide a trace of the phase response of a network at discrete 
frequency intervals. The group delay of the network is 
determined from the traces using the traditional method. 
Two points on the network analyzer's trace corresponding to 
the phase response of the network at two given frequencies 
are subtracted, then divided by the difference between the 
given frequencies. The result is an approximation of the 
derivative of the network's phase response as a function of 
frequency. The negative of the result is therefore an estimate 
of the group delay. The method can be repeated at discrete 
intervals along the phase response trace to generate a group 
delay trace. 

The traditional method of group delay estimation gener- 
ates two types of error from the true group delay that are 
dependent on the aperture width. First, noise or random 
measurement errors in the phase response measurements 
cause corresponding errors in the group delay estimates. The 
effect of the noise errors on the group delay estimates can be 
decreased by increasing the aperture width. For example, if 
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the aperture width is doubled and the noise variance of the 
phase response measurements is constant, the variance of the 
group delay estimates is halved. 

Larger aperture widths, however, have a smoothing effect 

5 on the group delay trace, resulting in a group delay trace 
which is distorted from the actual group delay. In places 
where the actual group delay is changing rapidly, for 
example, the group delay trace produced with a large 
aperture width is smoother than the actual group delay. This 

10 distortion from the actual group delay is a second type of 
error. So, large aperture widths are desirable for minimizing 
noise variance in the group delay trace, while small aperture 
widths are desirable for minimizing distortion of the trace; 
The selection of aperture size becomes a tradeoff between 

15 conflicting desired noise minimization characteristics to 
preserve the true shape of the group delay trace. 

The present invention provides a more accurate measure- 
ment of the group delay which has a decreased noise 
variance over a sequence of frequencies as compared to the 

20 traditional method. The present invention also provides a 
decreased trace distortion compared to the traditional 
method. In accordance with a preferred embodiment of the 
invention, the phase response of a network is measured at a 
sequence of stimulus signal frequencies. Using the 

25 unwrapped or cumulative phase response at a number of 
these frequencies within a frequency aperture, the phase 
response rate of change as a function of frequency is 
determined using linear regression techniques. Since the 
phase response rate of change as a function of frequency is 

30 the group delay, the result of the linear regression analysis 
yields a group delay estimate at a frequency in the center of 
the aperture. This group delay estimate determination is 
repeated at a sequence of frequencies to create a trace of the 
network's group delay across a range of frequencies. 

35 The method of the present invention is also applicable to 
estimating traces of measurements other than group delay 
that are the derivative of a measurement trace. For example; 
the method can be applied to estimate a trace of the 
instantaneous frequency of a signal over time from a trace of 

40 the signal's phase measured at discrete sample times. In a 
prior frequency estimation system, described in U.S. Pat; 
No. 4,983,906 to Hiller, a single frequency estimate is 
calculated using linear regression analysis techniques. The 
frequency estimate calculated by the Hiller system assumes 

45 that the frequency of the signal remains constant over time; 
However, it is often useful to monitor signals where the 
instantaneous frequency changes over time such as phase or 
frequency modulated signals or signals with drifting fre- 
quencies. With the present invention, it is possible to form 

50 a trace estimate of the instantaneous frequency over time. ; 
The foregoing and additional features and advantages of 
the invention will be more readily apparent from the fol- 
lowing detailed description, which precedes with reference 
to the accompanying drawings. 

55 DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of a network analyzer employ- 
ing a group delay estimator according to one embodiment of 
the present invention; 
6 0 FIG. 2 is a flow chart of a process according to one 
embodiment of the present invention; 

FIG. 3 is a flow chart of a computationally efficient 
process according to a second embodiment of the invention; 

FIG. 4 is a flow chart of a computationally efficient 
65 process for estimating a group delay trace from phase 
response sample deltas according to a third embodiment of 
the invention; 
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FIG. 5 is a graphical comparison of a group delay trace analysis is repeated to produce additional group delay esti- 

performed according to the present invention to a group mates at a sequence of other frequencies. The group delay 

delay trace performed according to the traditional method, estimates determined by the group delay processor 44 are 

both using the same aperture width and phase response data; indicated visually on a display 46. In some embodiments of 

and 5 the invention, the group delay estimates may also be stored 

™„ , . . - ji. ina data file or transferred to other computing devices for 

FIG. 6 is a graphical comparison of a group delay trace omer end X1SGS 

performed according to the present invention to a group w * 4l _ , , AA 

\ . e a \- ii . i More specifically, the group delay processor 44 applies 

delay trace performed conventionally on the same network. . r , ■ n/- i f .u c * * eA rr 

J r J the process shown in FIG. 2. In the first step 50, the group 

DETAILED DESCRIPTION OF THE 10 delay processor calculates a scale factor K used in the linear 

PREFERRED EMBODIMENT regression analysis formula described below. 

In step 52, the phase response processor sets a loop 

Referring to FIG. 1, a network analyzer 10 typically variable m equal to an initial value of N where 2N+1 is the 

comprises two input receivers, a response receiver 14 and a num ber of phase response samples within an aperture cen- 

reference receiver 16. When in use, the receivers 14, 16 are tered on lhe bin indexed by the loop variable m. (The 

connected to a network 18 under test. An input 22 of the aperture therefore includes N bins to each side of the bin 

response receiver 14 is connected to an output 24 of the indexed by tne loop varia51e m ) ^ fining steps 54, 56, 

network 18 to receive an output signal of the network. A 58 form a loop which is repe ated once for each value of the 

signal generator 26 produces a stimulus signal which drives variable m between N and L-N. At step 54 in the loop, the 

an input 28 of the network 18 and an input 30 of the group dday process0 r determines an estimate of the group 

reference receiver 16, The signal generator 26 is preferably delay d(m) at the frequency 0 f the bin indexed by the loop 

a component part of the network analyzer 10. However, an variable m. The loop variable m is incremented each itera- 

external signal generator can also be used as a source of the tion of the loop at step 56 Step 58 compares the loop 

stimulus signal. Typically, the stimulus signal is a sine variable m to the number L-N and repeats the loop until the 

waveform having a frequency which is swept through a loop variable exceeds L-N. Thus, the group delay processor 

predetermined range of test frequencies (i.e., a frequency wm dete rmine estimates of the group delay d(m) for the 

ramp signal). However, other forms of stimulus can be used mdices ioN to j^L-N. Calculating d(m) for only those 

such as noise, impulse or other arbitrary stimulus signals. indices betW een N and L-N ensures that there are an 

The output signal is produced by the network 18 in response adequate number of phase response samples to each side of 

to the stimulus signal. ^ the bin i ndexe d by the loop variable m to complete the 

Connected to the receivers 14, 16 are first-in, first-out frequency aperture. 

(FIFO) memory queues 36, 38. The receivers 14, 16 digitally In step 54, the group delay d(m) is calculated according 

sample the output and stimulus signals, producing streams t0 the following formula which is a standard linear regres- 

of digital samples of the signals which are placed in. the s j ori analysis formula in the form of the moving average 

memory queues 36, 38. 35 adap ted to the variables of phase response and frequency: ; 

The signal samples are read out of the memory queues 36, 

38 in FIFO order by a phase response processor 42. The (lj 

phase response processor 42 processes the digital samples of d(m) = K ^ (/• P[m +• fl) 

the stimulus and output signals in a well-known manner to i= ~ N 

produce phase response samples at uniform frequency inter- 40 

vals. The uniform frequency spaced, phase response samples where K (the scak f ac tor)-3/(N*(N+l)*(2*N+l)). 

are known as "bins," These phase response samples corre- xhe d(m) ^ equation (1) provides the slope of a 

spond to the unwrapped phase response of the network at phase ramp which is a least mean square error fit t0 the phasd 

various bins across a selected range of frequencies. The reS ponse samples within the aperture. However, the group 

phase response samples are "unwrapped" in the sense that 45 delay c f a ne twork is defined as the negative of the deriva-i 

the phase is not limited to the range 0 to 2n. Rather, the tive of the phase re sp 0 nse with respect to frequency! 

phase response samples correspond to the cumulative phase Therefore, the value given by the function d(m) in equation 

response of the network. (1) ^ actually the negat ive of an estimate of the network's 

The phase response samples are organized by the phase group delay. To properly estimate the group delay, the 

response processor 42 into an array, P[i], according to 50 function d(m) must be scaled by a negative one. Preferably; 

ascending order of stimulus signal frequency. The array P[i] this is accomplished by multiplying the scale factor K by a 

has indices i between zero and a maximum index number L. negative one so that K=-3/(N*(N+l)*(2*N+l)). The scale 

Each index i of the array corresponds to a uniformly spaced factor K may also be scaled by other factors to convert the 

frequency. The phase response sample or bin stored in the estimate to desired units of measure. Without further scaling; 

array at each index corresponds to the phase response of the 55 the function d(m) gives the group delay in radians or degrees 

network at the respective frequency. per bin depending on the units of the phase response 

The phase response samples are processed by the group samples. Thus, the scale factor K must be divided by the 

delay processor 44 to determine the group delay of the frequency spacing per bin in Hertz to determine group delay 

network 18 at a number of frequencies. In general, the group estimates in radians or degrees per Hertz. The scaling factor 

delay processor 44 applies linear regression analysis to 60 K may be scaled by other factors to obtain group delay 

determine a phase ramp having a least mean square error fit estimates in other units. 

to a set of the phase response samples. The set of phase The process described above is fairly computationally 

response samples used to determine the phase ramp on a intensive, averaging 2N addition operations and 2N+2 mul- 

selected frequency are those phase response samples within tiply operations per group delay estimate. In other words^ 

a frequency aperture centered on the selected frequency. The 65 each group delay estimate involves an "order of N". 

negative of the slope of this phase ramp corresponds to the calculation, a computational load that varies dependent on 

group delay of the network at the selected frequency. The the aperture width. Fortunately, regularities in the formula of 
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equation (1) can be exploited to lower the total number of d(N+l)) is determined in step 71 using equation (5). In the 

math operations required by the process where more than last step 72 of the initial portion of the process, a loop 

one group delay estimate as made. For example, a function variable m is set equal to N+2. 

f(m) equal to the difference between d(m) and d(m-l) can Next, the process has a loop portion comprising calcula- 

be defined as follows: 5 tion steps 74-80 and comparison step 82. In the first step 74 

of the loop portion, the value of the function g(m) is 

(j=m+N-i \ (2) calculated according to the equation (3) for the loop variable 

N (P[m + N]-P[m-N-l])- £ {P\J]) m. This value is used in step 76 to update the group delay 

Jsm ~ N ' delta accumulator. In step 76, the next group delay delta is 
10 calculated according to equation (4). The value g(m) from 

c ct \ • j i j n u * j ■ f step 74 is added to the value in the group delay delta 

The function f(m) is a group delay delta between adjacent ^ & c/ * / 

\ / o r j j accumulator to update the accumulator to f(m). This accu- 

1 A second function g(m) equal to the difference between ™ lalor value will be used in step 78 to update the group 

f(m) and f(m-l) can then be derived as follows: del T ay ^ ccuniu lator. . 

v 7 v 7 15 In step 78, the next group delay estimate is calculated 

«0»)«<K-N) (P[m+N]-P[m-N-2D + (KCN + i))(P[m-i-NhP[m-i + according to equation (5). In accordance with equation (5), 

N]) (3) the group delay delta f(m) for the frequency index m (the 

value in the group delay delta accumulator from step 76) is 

The group delay delta f(m) can also be expressed in terms of added to the value in the group delay accumulator to produce 

the function g(m) as follows: 2 0 the next group delay estimate d(m). This updates the group 

delay accumulator to the group delay for the frequency 

f[m)=Km-iyg(m) (4) mdex m> 

Jt _ iL j , ,* , ■%/ \ , , . In step 80, the loop variable m is incremented. Then in 

Further, the group delay estimate d(m) can be expressed in o^i * i_ i .\. . , . iU , T KT 

. „ r t . „*„ A~\t* t Aou* f^ii™,e- ste P 82, the loop variable is compared to the value L-N to 

terms or the group delay delta as follows: , ^ ,. lL r . , 4 t , . , . , . 

° 25 determine if the process is complete. If the loop variable m 

d(m)-d(m-i)+f{m) (5) ^ grater than L-N, the process is complete. Otherwise, the 

loop portion of the process is repeated using the incremented 

The functions f(m) and g(m) can be used to implement a loop variable m. 

more computationally efficient process for estimating a As each group delay estimate is calculated in step 78, it 

plurality of group delays at a sequence of frequencies. Using 30 is stored in a group delay trace preferably in the form of an 

these equations, the group delay estimates can be made by array. Thus, the process requires that data storage for both an 

a pair of accumulators. A first or group delay accumulator array of the phase response samples P[i] and an array of 

stores the current group delay estimate d(m). A second or values of the function d(m) for m=N to m=L-N are main^ 

group delay delta accumulator stores the group delay delta tained. However, the only values of the functions f(m) and 

f(m). After determining the initial group delay estimate and 35 g(m) that need be maintained are those for the current 

group delay delta, subsequent group delay estimates can be calculation of d(m) within the loop portion (steps 74-82) of 

determined by, first, calculating the difference function g(m) the process. By comparison, in the process shown in FIG. 2, 

between the current and next group delay deltas according to two arrays for storing the phase response trace and the group 

equation (3). Then, updating the group delay delta in the delay trace are also maintained. The computationally inten- 

second accumulator by adding the calculated value of g(m) 40 sive process shown in FIG. 3, therefore, requires minimal 

for the next frequency. Finally, the group delay estimate is additional data storage. 

updated in the first accumulator by adding the group delay A further modification to the above processes can be madd 

delta from the second accumulator. In this manner, only five to avoid problems of numerical representation. In the pro-; 

addition or subtraction operations and two multiplication cesses described above, the unwrapped phase response is 

operations are required to determine each subsequent group 45 used. If the phase response of a network varies over many 

delay estimate. With a fixed number of operations per cycles in the tested frequency range, a large number of digits 

subsequent group delay estimate, each subsequent estimate will be required to represent the untrapped phase response 

is an order of one calculation and invariant of changes in the with sufficient precision. For example, if four hundred phase 

aperture width. Thus, for large aperture widths, the compu- response samples are obtained and the difference in phase 

tationally efficient method offers significant performance 50 response between adjacent samples is nearly 180 degrees 

improvement. (the maximum difference in phase response between adja- 

F1G. 3 illustrates in more detail this computationally cent samples that can be detected unambiguously), then 

efficient form of the process employing the f(m) and g(m) nearly 200 cycles of phase response will have been accu- 

f unctions. The initial portion of the process comprises steps mulated between the first and last phase response samples: 

66-72. The initial step 66 in the process is to calculate the 55 About 9 bits of mantissa are required to represent this 

scale factor K used by the functions d(m), f(m), and g(m). amount of unwrapped phase response with sufficient precis 

Preferably, the values of K*N and K(N+1) are also calcu- sion. 

lated since they are constants used by the function g(m) for To avoid numerical representation problems caused by th4 

each group delay estimate. In the second step 68, the group use of the unwrapped phase response, the processes 

delay processor determines the group delay estimate for the 60 described above can be modified so that the phase response 

frequency index N according to equation (1), i.e. d(N) is samples are represented in terms of the difference in phase 

determined. This determines the initial value of the group response from a previous adjacent sample (phase response 

delay accumulator. In step 70, the value of f(N+l) is delta) instead of the unwrapped phase response. The linear 

determined according to equation (2). The value f(N+l) is regression formula of equation (1) must be modified to 

the initial value for the group delay delta accumulator. Using 65 account for this change. With the phase response delta 

the values for d(N) and f(N+l) determined in steps 68, 70, samples organized into an array D[i] such that, for i«l to 

the group delay estimate for the frequency index N+l (i.e., i-L: 
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DlO-Pl/J-Pli-i] (6) ment of the group delay of a network over a range of 

frequencies between 1.6 and 1.8 MHz using the traditional 
the formula is as follows: method of the pfior ^ A tface n is a measurement 

of the group delay of the same network over the same 
j/ v K V /a, i ^ m ™ ^ 5 frequency range applying the present invention. An aperture 

2 containing 64 phase response samples was used for each 

group delay estimate in the traces 90, 92. The trace 90 using 
the traditional method has a significantly greater variance 
The equation (7) can be used to compute each group delay tnan tne trace 92 us ing the present invention. Theoretically, 
estimate directly in a process identical to that shown in FIG. 10 f or a i arge number n of samples in the aperture, the decrease 
2. A computationally efficient process similar to the process i n variance provided by the present invention should be 
shown in FIG. 3 can also be derived from the equation (7). approximately n/6. In this example where there are 64 
However, one additional accumulator step is required. A samples in the aperture, the variance of trace 92 should be 
group delay delta function r(m) equal to the difference about n tj mes i ess t h an the variance of trace 90. 
between d(m) and d(m-l) is derived as follows: 15 Ifl FIG 6) a group dday trace m of a network over a 

range of frequencies from 5 MHz to 125 MHz performed 
( ^ ] W according to the present invention is compared to a group 

*m) = K-y^ (i-D(m + 0)j y trace m Qf ^ same netwQrk per f ormed the 

traditional method of the prior art. The group delay trace 112 
20 of the present invention has a peak group delay 116 higher 
A function s(m) equal to the difference between r(m) and than a peak group delay 118 of the trace 114. The lower peak 
r(m-l) is derived as follows: group delay 118 is due to the smoothing effect of the 

traditional method. The present invention produces less 

4<=«^i } (9) distortion than the traditional method, resulting in the group 

/v ■ (D{m = N) + D[m - 1 - N)) - £ (D(m + 2 5 delay trace 112 that is closer to the actual group delay of the 

,= ' N network. 

The method of the present invention is also applicable to 
A function t(m) equal to the difference between s(m) and estimating traces of measurements other than group delay 
s(m-l) is derived as follows: that are the derivative of a measurement trace. For example; 

30 the method can be applied to estimate a trace of the 
^S^riS^ " D(m " N " 2))+K(N+1) ^ W " 1 - N >- instantaneous frequency of a signal over time from a trace of 

the signal's phase measured at discrete sample times. Either 
The function s(m) can then be expressed in terms of the the form of the method illustrated in FIG. 2 or the form 
function t(m) as follows: illustrated in FIG. 3 can be used. The only required modi- 

35 fication to the methods is to rescale the constant multipliers 
sim)=s(m-i)+t(m) (11) j n me e q Ua tions so that frequency, which is the rate of 

The function n» expressed in terms of the function s(m) is: chan S e of P hase versus time ' * calculated. Further, if the 

phase trace is in the form of phase delta samples (i.e., each 
r(m)-/<m-i)+j(m) (12) sample in the trace is the difference in phase from the 

. ,/ v , . r.x. c 40 previous sample time), the methods for performing an 

Finally, the function d(m) expressed in terms of the function r . c JU ♦ u *u *u j *n » ♦ a • '■ 

( . / \ \ / f estimate from a delta trace such as the method illustrated in 

r{m) is. ^ can ^ usec j Again, the constant multipliers in the 

d(m)=d(m-i)+r(m) (13) equations must be rescaled. 

„ T . , „ A . - . / v , x * Having described and illustrated the principles of my 

With reference to FIG. 4, the functions r(m), s(m), and 45 ^ reference tQ a ferred embodiment ^root, 

t(m) are used in a computationally efficient process for ft ^ be t ^ ^ invention can be modified m 

estimating the group delay trace from the phase response arrangemem and detail without departing from such pru> 
delta which xs sumhr to the ; computationally efficient pro- d , es For h ^ fa in FJG 1 ^ 

cess shown in FIG. 3 In the process shown m FIG. 4, e base ^ an(J del lhe 

however an additional accumu ator is required, i.e the 50 invention ^ also be irap i emented ^ a single processor 
functions d(m), «m), and s(m) all operate as accumulators suiuble fof executing the processes carried out by toth of 

updated in turn according to the equations (llWlJ). More _ . 

r ii • . v , , 1ftft ( ? , v p , , the separate processors, 
specifically, in an initial step 100, the values of the constants _ , , .„ , . , , . . , ,. j 

K, K*N, and K(N + 1) are calculated. In step 102, the value . Further > * ltho "S h ll,u k strated Wlth the 1 " ve " tl0n embod , ied 

of d(N) is determined according to equation 7. In step 103, 55 in u a netw ° rk ^ mT ^ lhe inv f ntl0n can also be em P Io y ed in 

the value of r(N + l) is determined according to equation (8) ° ther ^cations not involving a network analyzer tp 

then d(N + l) is determined from r(N + l) using equation (13). determme the ^ of a network ^ en P r0Vlded Wlth 

In step 104, the value of s(N + 2) is determined according to measurements of the phase response of the network, 

equation (9) then r(N+2) and d(N+2) are determined from T° e ^rm of the methods can also be modified without 

s(N+2) using equations (12) and (13), respectively. The 60 departing from the principles of the invent,on. Among other 

values of d(m) for m«N+3 to m=L-N are determined within modifications in the method illustrated in FIG. 3, steps 74 

a loop comprising steps 105-108. In step 106 in the loop, the a °d 76 of calculating g(m) and f(m) need not be performed 

value of t(m) is determined according to equation (10) then separately. The calculation of g(m) can be done implicitly in 

the values of s(m), r(m), and d(m) are determined in order a direct update of f(m). 

according to equations (11)-(13), respectively. 65 Accordingly, I claim as my invention all such modifica- 

FIG. 5 shows an example of the results obtainable through tions as may come within the scope and spirit of the 

use of the present invention. A first trace 90 is a measure- following claims and equivalents thereto: 
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I claim: 

1. In an electronic measurement instrument, a method for 
estimating a trace of second measurements related to the 
derivative of a trace of first measurements of an electrical 
property of a network, the method comprising: 

obtaining a plurality of uniformly spaced samples of a 
first measurement to form the first trace; and 

performing a least mean squares fit to samples in the first 
trace within each of a plurality of apertures to deter- 
mine a plurality of estimates of a second measurement 
for forming the second trace. 

2. The method of claim 1 further comprising multiplying 
the second trace by a constant multiplier to form a trace of 
estimates of the second measurement. 

3. The method of claim 2 wherein the first measurements 
are the phase response of a network sampled at uniformly 
spaced frequencies and the second measurements are the 
group delay of the network at various frequencies. 

4. The method of claim 1 wherein the first measurements 
are the phase of a signal sampled at uniformly spaced sample 
times and the second measurements are the instantaneous 
frequency of the signal at various times. 

5. The method of claim 1 comprising: 

determining a first estimate in the second trace by per- 
forming a least mean square error fit to samples of the 
first trace within a first aperture; 

determining each subsequent estimate in the second trace 
by determining a value which, when added to a previ- 
ous estimate adjacent the subsequent estimate, provides 
a least mean square error fit to samples of the first trace 
within a subsequent aperture of the plurality of 
apertures, and adding the value to the previous estimate 
to form the subsequent estimate. 

6. An electronic measurement instrument comprising: 

a reference receiver connectable to an input of a network 
to receive a stimulus signal; 

a response receiver connectable to an output of the 
network to receive a network output signal produced by 
the network in response to the stimulus signal; 

means for measuring the phase response of the network at 
uniform intervals of frequency; and 

means for performing a least mean squared error analysis 
on a number of phase response measurements made at 
frequencies within each of a plurality of apertures to 
estimate a group delay trace of the network at a 
plurality of selected frequencies. 

7. The instrument of claim 6 wherein the means for 
measuring the phase response of the network is a phase 
response processor for determining an array of phase 
response data measured at uniform frequency intervals from 
samples of the stimulus signal and the network output signal. 

8. The instrument of claim 7 wherein the means for 
performing the least mean square error analysis is a group 
delay processor for performing the least mean squared error 
analysis on subsets of the array to yield estimates of the 
group delay of the network at a plurality of selected fre- 
quencies. 

9. The instrument of claim 6 wherein the means for 
measuring the phase response and the means for performing 
the least mean square error analysis are embodied in a single 
processor. 

10. The instrument of claim 6 wherein the phase response 
measurements are measurements of the cumulative phase 
response of the network. 

11. The instrument of claim 6 wherein the phase response 
measurements are measurements of the difference in phase 
response of the network between two adjacent frequencies. 
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12. The method of claim 3 comprising: 

determining a first estimate in the second trace by per- 
forming a least mean square error fit to samples of the 
first trace within a first aperture; 

determining each subsequent estimate in the second trace 
by determining a value which, when added to a previ- 
ous estimate adjacent the subsequent estimate, provides 
a least mean square error fit to samples of the first trace 
within a subsequent aperture of the plurality of 
apertures, and adding the value to the previous estimate 
to form the subsequent estimate. 

13. The method of claim 3 further comprising multiplying 
the second trace by a constant multiplier to form a trace of 
estimates of the second measurement. 

14. The method of claim 4 comprising: 

determining a first estimate in the second trace by per- 
forming a least mean square error fit to samples of the 
first trace within a first aperture; 

determining each subsequent estimate in the second trace 
by determining a value which, when added to a previ^ 
ous estimate adjacent the subsequent estimate, provides 
a least mean square error fit to samples of the first trace 
within a subsequent aperture of the plurality of 
apertures, and adding the value to the previous estimate 
to form the subsequent estimate. 

15. The method of claim 4 further comprising multiplying 
the second trace by a constant multiplier to form a trace of 
estimates of the second measurement. 

16. In an electronic measurement instrument, a method 
for estimating a trace of second measurements related to the 
derivative of a trace of first measurements of an electrical 
property of a network, the method comprising: 

obtaining a plurality of uniformly spaced samples of a 
first measurement to form the first trace; and 

performing a least mean squares fit of a linear function to 
samples in the first trace within each of a plurality of 
apertures, wherein at least two of the apertures have 
more than one of the samples in common; and 

determining a plurality of estimates of a second measure 
ment centered in the apertures to form the second trace 

17. The method of claim 16 further comprising multiply 
ing the second trace by a constant multiplier to form a trace 
of estimates of the second measurement. 

18. The method of claim 16 wherein the first measurer 
ments are the phase response of a network sampled at 
uniformly spaced frequencies and the second measurements 
are the group delay of the network at various frequencies.; 

19. The method of claim 16 wherein the first measurer 
ments are the phase of a signal sampled at uniformly spaced 
times and the second measurements are the instantaneous 
frequency of the signal at various times. 

20. The method of claim 16 comprising: 
determining a first estimate in the second trace by per^ 

forming a least mean square error fit to samples of the 
first trace within a first aperture; 
determining each subsequent estimate in the second trace 
by determining a value which, when added to a previa 
ous estimate adjacent the subsequent estimate, provides 
a least mean square error fit to samples of the first trace 
within a subsequent aperture of the plurality of 
apertures, and adding the value to the previous estimate 
to form the subsequent estimate. 



07/09/2004, EAST version: 1.4.1 



